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ABSTRACT 

*  ^By  means  of  nonlocal  elasticity  theory,  the  Griffith  crack  problem 
Is  solved  to  determine  the  state  of  stress  In  a  plate  weakened  by  a  line 
crack.  The  maximum  stress  Is  found  to  be  finite  and  occurs  away  from  the 

*  crack  tip.  Cohesive  stress  is  estimated  and  compared  with  the  calculations 
based  on  atomic  lattice  theory. 

/V 

1.  INTRODUCTION 

In  several  previous  papers,  we  have  Investigated  the  state  of  stress 
In  plates  containing  sharp  line  cracks  [1,  2,  3]  and  in  Volterra  dislocations 

* 

[4,5].  These  and  our  other  similar  work  Indicated  that  the  nonlocal  theory 
eliminates  the  stress  singularity  predicted  by  the  classical  elasticity  theory. 
These  results  enabled  us  to  introduce  a  natural  fracture  criterion  based  on 

t  p 

maximum  stress  hypothesis  and  calculate  the  theoretical  cohesive  strength  of 
various  crystalline  materials. 

In  reference  [1],  the  maximum  stress  was  predicted  to  occur  at  the 
crack  tip  while  the  solution  of  screw  dislocation  problems  showed  that  the 
maximum  stress  occurs  slightly  away  from  the  edge  of  the  dislocation.  Accord¬ 
ing  to  a  lattice  model  of  crack  used  by  Elliot  [6],  the  maximum  stress  occurs 

I 

slightly  away  from  the  crack  tip.  This  situation  was  also  confirmed 
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In  the  recent  computer  simulations  of  Gohar  [7].  To  obtain  such  fine  results 
.9  near  the  crack  tip,  computer  calculations  must  Include  discretization  of  a 

dimension  less  than  one  atomic  distance  and  must  also  overcome  the  major  dif¬ 
ficulty  related  to  the  round-off  errors.  In  [1],  such  fine  resolution  was 
^  not  Introduced  and  It  prevented  the  prediction  of  the  exact  location  of  the 

maximum  stress.  Moreover,  In  [1],  two-dimensional  generalization  of  a  one- 
dimensional  nonlocal  kernel  was  used. 

>  In  classical  elasticity  theory,  a  mixed-boundary  value  problem, 

prescribing  the  tractions  on  one  part  of  the  boundary  surface  and  the  dis¬ 
placement  fields  on  the  remainder.  Is  well-posed.  In  nonlocal  elasticity, 

^  this  Is  not  the  case  because  of  the  fact  that  the  prescription  of  the  stress 

field  at  one  point  requires  the  knowledge  of  the  strain  field  in  the  neigh¬ 
borhood  of  the  point.  Since  at  the  crack  tip  this  region  extends  into  the 
^  body.  It  Is  necessary  to  verify  the  self-consistency  of  the  boundary  condi¬ 

tions  In  this  overlapping  region.  Without  such  consideration,  the  boundary 
conditions  within  two  atomic  distances  around  the  crack  tip  Is  satisfied  in 
t  a  non-uniform  approximate  manner  (Atkinson  [8]). 

The  raison  d'etre  of  the  present  work  takes  its  origin  from  the 
foregoing  considerations.  In  Section  2,  we  introduce  a  two-dimensional  kernel 
*  by  matching  the  dispersion  curve  of  a  square  lattice  to  that  of  nonlocal  plane 

waves.  In  Section  3,  a  set  of  self-consistent  boundary  conditions  are  rederived. 
Section  4  deals  with  the  solution  of  the  Griffith  crack  problem.  Analytical 
I  expressions  are  given  for  the  stress  field  along  the  crack  line.  The  maximum 

stress  and  Its  location  are  calculated.  In  Section  5  we  compare  these  results 
with  those  known  In  the  atomic  lattice  theory  and  computer  simulation  studies. 

1  Quantitative  agreements  are  gratifying. 
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2.  TWO-DIMENSIONAL  KERNEL 

From  [9],  for  a  two-dimensional  square  lattice,  we  have  the  dispersion 
relations 


(2.1)  or  *  $q{2  -  cos(Ea)  -  cos(na)  +  2  [1  -cos(?a)  cos(na)]) 


(2.2)  a)  *  eQ[2  -  cos(?a)  -  cos(na)] 

where  £  and  n  are  Fourier  variables  corresponding  to  x  and  y  directions, 
a  Is  atomic  lattice  parameter  and 


(2.3)  eQ  =  f^/a)2  , 


C]  =  [(X  +  2y)/p] 


Here,  c-j  is  the  phase  velocity  of  the  irrotatlonal  waves,  (X,y)  are 
Lame  constants  and  p  Is  the  mass  density. 

In  nonlocal  elasticity  for  isotropic  solids,  the  stress  constitutive 
equations  are  given  by  [10], 


(2.4)  t. 


■  f  tX'dx'-xl)  err(x-)«kl,  +  ai-dx'-xDe^fx-ndaCx') 


where  the  integral  Is  over  the  two-dimensional  plane  region  R(x',y')  and 
ekt  the  strain  tensor  which  Is  related  to  the  displacement  vector  u^  by 

(2.5)  «kt-  I(uk.t  +  Uil.k) 

Here,  and  throughout  an  Index  following  a  comma  represent  gradient. 
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Uk.£  *  3uk/3xJt 

In  this  section  we  determine  the  kernels  X1  and  y'  by  comparing 
(2.1)  and  (2.2)  with  the  dispersion  of  plane  waves  In  nonlocal  elasticity. 
Substituting  (2.4)  and  (2.5)  Into  the  equations  of  motion 


> 


» 


<z-6>  ‘kt.k  -  0“t  *  0 

and  taking  the  two-dimensional  Fourier  transform  of  the  ensuing  field 
equations,  we  obtain 

[(A,+2{i,)S2  +  y'n2  -  pw2)]0  +  (X*  +  y')  £n  v  =  0 

(2.7) 

(A'  +  u 'fen  U  +  [(A*  +  2y')n2  +  jj 'C2  -  pw2]  v  *  0 


> 


» 


I 
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where  a  superposed  bar  denotes  Fourier  transform 


(2.8)  u(C,n) 


00  00 

2ir 

•  00  *00 


u(x,y)  exp[1(cx  +  ny)]  dxdy 


For  a  non-trlvlal  solution  of  (2.7)  to  exist,  the  determinant  of 
the  coefficients  of  u  and  v  In  (2.7)  must  vanish,  leading  to  dispersion 
relations. 


(2.9) 

(2.10) 


“2/c2  -  (C2  +  n2)(X’+  2u* ) /(X  +  2y) 
w2/c2  *  U2  +  n2)  S'  /y 


» 


Comparing  (2.9)  with  (2.1)  and  (2.10)  with  (2.2),  we  obtain  expressions  for 
X’  and  y'  . 

Eventually,  our  main  Interest  Is  In  the  nature  of  the  singularity 
of  X'  and  y1  .  An  examination  X1  and  y1  so  obtained  showed  that 
they  both  possess  logarithmic  singularities.  In  fact.  Inverse  transform 
of  y*  yields,  for  the  case  of  Poisson  material  (X  »  y): 

(2.11)  y'(x»y)  -  2a-2  in  {[(x-a)2  +  y2][(x+a)2  +  y2][x2+ (y-a)2] 

[x2  +  (y+a)2]  (x2  +  y2)“4)  ** 

It  Is  clear  that  y'  is  not  an  Isotropic  kernel.  However,  noting  that  It 
has  a  logarithmic  singularity,  we  may  construct  an  isotropic  kernel  with  the 
same  singularity  for  both  X'  and  y'  .  The  following  kernel  satisfies 
the  basic  requirements  of  a  two-dimensional  kernel. 

(2.12)  o(|x3-x|)  *  ^  B2K0{e[(x'-x)2  +  (y’-y)2]*} 

where  KQ  is  the  zeroth  order  modified  Bessel  function,  0  Is  a 
parameter,  and 

a  =  X'/X  *  y'/y  . 

He  note  the  following  basic  properties  of  o: 

a)  It  has  logarithmic  singularity  at  the  origin 

b)  It  vanishes  as  Ix'-xl  ® 

m*  m*m* 

c)  It's  Integral  over  the  plane  Is  unity,  l.e. 

f  r 


(2.13) 


o(|x'-xj  da(x') 
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d)  As  B  ♦  •  •  It  reverts  to  a  Dirac  delta  measure,  so  that 
(2.4)  gives  classical  Hookes  law. 

Also,  this  kernel  possesses  an  Interesting  property  In  that  it  satisfies 
the  following  differential  equation 

(2.14)  +  r  If  "  6?)  Ko(6r)  *  s2fi(r) 

where  6(r)  Is  the  Dirac  delta  function. 

This  property  Is  very  convenient  for  the  reduction  of  the  Integro- 
dlfferentlal  field  equations  of  the  theory  of  nonlocal  elasticity. 

To  determine  @  ,  we  match  the  dispersion  curve  (2.2)  at  the  end 
of  Brlllouln  zone.  This  was  done  for  two  extreme  cases. 


(2.15) 

1/6 

«  0.22  a 

(C 

*  ir/2a  , 

n  *  w/2a) 

(2.16) 

1/6 

*  0.31  a 

(5 

=  ir/2a  , 

n  ■  0) 

The  maximum  error  In  these  cases  are  less  than  1%  and  0.5%,  respectively. 

By  Integrating  o  over  y'  from  -»  to  «  we  can  also  derive 
a  one-dimensional  kernel  which  Is  of  the  form 

(2.17)  o(|x'-x|)  «  |  B  exp[-B  |x'-x|] 

It  must  be  noted  that  none  of  these  kernels  possess  finite  support,  there' 
fore  special  attention  should  be  paid  to  the  asymptotic  behavior  of  the 
solutions  Involving  Fourier  Integrals. 


3.  BOUNDARY  CONDITIONS 


A  plate.  In  the  (x,y)  plane,  containing  a  line  crack  of  length 
2 1  along  the  axis  of  the  rectangular  coordinates.  Is  subjected  to  a  uniform 
tension  tQ  at  y  =  »,  in  the  y-dlrectlon  (Figure  1).  This  problem, 
which  constitutes  the  foundation  of  contemporary  fracture  mechanics.  Is 
known  as  the  Griffith  crack  problem,  first  treated  by  Griffith  [11].  In 
the  classical  elasticity  treatment, appropriate  boundary  conditions  are: 

*  0.  y.o,  V« 

(3.1)  ‘yy  •  -  V  y  ■  0  •  1*1  *  1 

v  *  0  ,  y  *  0 ,  |x|  >_  £ 

where  t££  represents  the  classical  stress  field.  The  Griffith  crack  problem 
requires  the  addition  of  a  uniform  stress  field  t*L  *  tQ  to  the  solution 
obtained  under  the  boundary  conditions  (3.1).  These  boundary  conditions  do 
not  carry  over  to  the  nonlocal  elasticity  in  a  strict  sense.  This  is 
because  at  the  crack  tips  (x  *  ±  l  ,  y*0)  both  the  traction  and  the 
displacement  fields  cannot  be  prescribed  simultaneously.  Because  of  the 
Integral  form  (2.4)  of  the  stress  constitutive  equations  prescription  of  the 
stress  field  at  the  crack  tip  requires  that  strain  field  must  be  known  In  a 
region  around  the  tips.  Due  to  the  short  range  of  Interatomic  force,  this 
range  Is  of  the  order  of  few  atomic  distances.  Nevertheless,  for  proper 
continuum  formulation.  It  Is  necessary  to  state  the  boundary  conditions  in 
a  self-consistent  form.  To  this  end  we  observe  that 


I  ,in.  nmmwmmu  ..ui  win" 

l 

(a)  In  the  limit  as  8-*«>  the  nonlocal  stress  field  must 
approach  to  the  local  stress  field. 

*  (b)  The  stress  field  must  be  continuous  In  B  for  all 

x  and  y  and  possess  continuous  partial  derivatives 
with  respect  to  6  for  all  6  . 
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The  stress  constitutive  equations  (2.4)  of  the  nonlocal  theory  may 
be  expressed  as 

(3-2)  lkt  ■  jj  '“(2’) 

R 

where  a Is  the  local  stress  tensor  given  by  the  Hookes  Law: 

°k*  *  Xerr6kt  +  2liekt 
If  we  express  In  the  asymptotic  form 

(3.3)  a(x,0)  *  f n(x)  +  l  f  (x)  3"n 

~  *»  ~  u  ~  ~n  ~  . 


Using  (3.3)  In  (3.2)  and  applying  the  nonlocal  smoothness  conditions 

t 

(3.4)  11m  8^+n  £— •  (8n"^  t  (x,B)}  1  0  ;  n>0 

30"  ~  ~ 

•  we  obtain  restrictions  on  the  functions  f„(x)  . 

-n  +* 

For  the  Griffith's  problem,  this  leads  to 


(3.5) 


f°J-  '  '  *0  •  I 

f0  (x)  -0  1*1  ‘  ‘ 

UXj  ~ 
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and  all  other  f„(x)  vanish  Inside  of  the  crack  line. 

Boundary  conditions  of  the  nonlocal  theory,  appropriate  to  the 
Griffith  problem,  are  therefore  given  by 


’  °xy(x,0)  *  0  • 

y  *  0  ,  V  x 

(3.6) 

t‘y(*.0) 

“  Oyy(X.  0)  =  *  tQ  , 

y  =  o ,  |  x  |  <  t 

V 

=  0 

y *  o  |x|  _>  i 

Note  that  these  boundary  conditions  are  Identical  to  those  of  the  classical 
theory  (3.1).  This  Is  due  to  the  fact  that  the  applied  tension  tQ  ,  along 
the  crack,  is  constant.  For  a  non-uniform  loading,  tQ  In  (3.6)  Is  replaced 
by  a  more  complicated  function  and  the  functions  f_  may  not  vanish.  As  a 

Mil 

result,  boundary  conditions  will  be  complicated  by  the  requirement  that 
t£  =0  is  no  more  valid  along  the  entire  x-axis. 

Xjr 
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4.  THE  SOLUTION 

Following  Eringen,  et  al  [1],  for  the  field  equations  of  the  two- 
dimensional  nonlocal  elasticity,  we  have 


ff  a(|x'-x|)[(X+y)(u,  .  .  +  y^2u']dx'dy’ 

R  x  x  x  y 

l 

-2 y  /  a(|x’-x|){e'  (x\0)]dx*  -  0 

-i  y* 

//o(|x*-x|)[(X+y)(u'  ,  ,  +▼!_,_,)  +  yvVldx'dy’  -  0 
R  x  y  7  7 
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where 


(4.2)  [e^x’.O)]  '  *yx(x,,0+>  “  V(*,,0  ) 


For  the  kernel  we  use  (2.12).  Applying  the  operator  In  (2.14) ,  using  the 
Fourier  transform  along  the  x-directlon 


(4.3) 


f(k.y) 


f(x,y)  exp(ikx)  dx 


-  oo 

we  obtain  the  general  solution  of  (4.1) 


u  -  (2*)'*  /  i  [jk|  A(k)  +  (Jkjy  -  TJ~)B(k)  lerpHkly  -  ikx)dk 

(4.4)  -  L  **  J 

v  -  (2v)  **  f  jA(k)  +  yB(k)  Jexp(-|kjy  -  lkx)dk 

The  boundary  condition  (3.6)  leads  to  a  pair  of  dual  integral 
equations  for  A(k)  whose  solution  determines  A(k)  and  B(k)  : 


(4.5)  A(k)  *  A0J}(k)/k 
B(k)  *  Bq  |k|  A(k) 

where 

Aq  =  (2/tt)1/2  [2p  (\+v)/l2  tQ(X+  2p)] 

(4.6)  Bq  *  (A+y)/(X  +  2m) 

k  B  kt 


I 


t 


Upon  substituting  (4.5)  Into  (4.4)  and  (3.2),  we  obtain  displacement 
and  stress  fields.  Relevant  to  the  present  work  Is  the  hoop  stress  along 
the  crack  line,  which  Is  given  by: 


OCX 


1 

(4.7) 

tyy(5,0)/tQ  =  -  cx(ek)  J^tc)  cos  (<C)  die 

» 

where 

(4.8) 

. 

0 

S  =  x/t  ,  e  *  1/Bit 

(4.9) 

a(e<)  =  [(1  +  e2ic2)^  +  ck]~^  +  ek[(1  +  e2ic2)^  +  ek]-^ 

'  f 

r 

\ 

\ 

I 

a(eic) 

given  by  (4.9)  is  much  simpler  than  (5.11)  obtained  in  [1], 

f 

;  * 

The  integral  in  (4.7)  was  evaluated  by  use  of  the  asymptotic  expansion 

; 

of  J.|  (k)  : 

(4.10) 

Jj(k)  n.  k-V2  cos(k  _  3tt/4)  +  0(k  ^)  ;  |k |  -*-  <» 

t 

The  result  to  the  order  one  (e°),  is 

I 

i 

S 

(4.11) 

1 i2  \V2 

tyy(£,0)/tQ  *  (tt/2)  {-g—  [I1/4  *3/4  +  I3/4  K|/4] 

1 

*  | 

*1/2 

+  2e  ^1/4  K3/4  '  1 5/ 4  K7/4  +  X3/4  Kl/4  "  X7/4  K5/4 

i 

where 

(4.12) 

X  *  £  -  1 

»  and  Iv,  are  order  modified  Bessel  functions  and  they  are  evaluated 
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At  the  crack  tip  C*1  and  we  obtain 


(4.13) 


t  (l.O/tfl^  0.57446  e' 


We  notice  however,  that  the  hoop  stress  Is  not  maximum  at  the  crack 
tip.  In  Table  1  we  list  the  hoop  stress  In  the  vicinity  of  the  crack  tip 
as  a  function  of  A/2e  .  From  this  table  It  Is  clear  that 


(4.14) 


Vy«*  *  '-'2  V'-0! 


and  this  occurs  at 

(4.15)  5*1+  0.50  e 

Using  the  values  of  e  given  (2.15)  and  (2.16),  we  obtain 


» 

(4.16) 

V  max^o  *  0-8'  • 

5  *  1  +  0.1  a/l 

(4.17) 

lyy  max's  ’  0.97  (Zt/.)1'2  . 

5  %  1  +  0.15  i/l 

5.  COMPARISON  WITH  ATOMIC  LATTICE  THEORIES 

The  hoop  stress  along  the  crack  line,  given  by  (4.11),  is  plotted 
In  Figure  2  against  5  *  x/£  for  e  *  0.22  a/£  .  On  this  figure,  the 
hoop  stress  obtained  by  Elliot  [6],  and  Errngen  et  al.  [1]  are  also  shown. 
The  maximum  stress  and  Its  location  and  the  ratio  of  the  maximum  stress  to 
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the  stress  at  the  crack  tip  are  compared  In  Table  2.  In  Elliot's 
I  calculations,  21/ &  -  1000  and  the  Poisson's  ration  v-0.25.  We  have  there¬ 

fore  used  the  same  numbers  In  (4.16)  and  (4.17)  and  the  quantitative 
agreement  Is  very  good. 

I  Gohar's  [7]  recent  computer  simulations  on  the  crack  problem  In  a 

two-dimensional  triangular  lattice  Indicates  that  the  hoop  stress  profiles 
are  of  the  same  character  as  those  of  reference  [1]  and  the  present.  He 
^  found  that  *s  9enera11y  larger  than  nonlocal  continuum 

results  while  his  stress  profiles  were  somewhat  broader.  Strict  comparison 
with  his  results  may  not  be  meaningful  because  of  the  Inherent  nonlinearity 
(  In  the  force  potentials  used  In  his  model.  However,  In  some  cases,  the 

maximum  stress  and  the  stress  profile  were  found  to  be  In  good 
quantitative  agreement  with  the  nonlocal  model. 

^  Cohesive  stress  can  be  calculated  by  use  of  (4.16)  or  (4.17)  In  the 

same  manner  as  In  [1  -  5],  If  one  uses  (4.16),  one  obtains  10%  higher 
cohesive  stress  and  if  (4.17)  is  used,  the  cohesive  stress  becomes  30% 

^  higher  than  those  given  in  reference  [1], 

The  nonlocal  fracture  criteria  Introduced  by  Erlngen  [12]  can  be 
used  in  the  same  manner  by  setting  t|flax  =  Cohesive  * 

?<  The  present  work  modifies  the  boundary  conditions  used  previously 

in  [1]  to  eliminate  overdetermination  arising  from  the  overlapping  Influence 
region  of  the  stress  and  displacement  fields  on  the  boundary.  Introduction 
I  of  a  new  two-dimensional  kernel  brings  the  nonlocal  model  closer  to  two- 

dimensional  lattice  models.  Nevertheless,  numerical  results  are  very  close 
to  those  In  [1],  probably  because  the  Influence  of  the  boundary  "overlapping 
I  region"  Is  extremely  small,  not  affecting  the  outcome  appreciably.  The  present 

approach  enables  us  to  state  that  the  maximum  stress  does  not  occur  at  the  crack  tip. 
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TABLE  1 

HOOP  STRESS  HEAR  THE  CRACK  TIP 


A/2e 

tvy(e,o)/t, 

0.0 

1. 

0.1 

1.09 

0.2 

1.10 

0.25 

1.12 

0.3 

1.10 

1.0 

0.89 

1.5 

0.75 

2.0 

0.66 

TABLE  2 

HOOP  STRESS  ALONG  THE  CRACK  LINE 
_ W^O _ ^max _ ^frax^tlp 


Elliot  [6] 

27.62 

1  +  0.2  (a/£) 

1.16 

Nonlocal  Ref.  [1] 

23.08 

— 

— 

Nonlocal  |  1/Ba  *  0.31 

25.61 

1  +0.15  (a/A) 

1.12 

Present  i  1/Ba  ■  0.22 

30.67 

1  +  0.1  (a/A) 

1.12 
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